CERN-PH-TH/2007-075 
IFT-UAM/CSIC-07-19 



Sterile neutrino oscillations after first MiniBooNE results 

Michele Maltoni^'0 and Thomas Schwetz^'0 

^ Departamento de Fisica Tedrica & Instituto de Fisica Teorica, 
Facultad de Ciencias C-XI, Universidad Autonoma de Madrid, 
Cantohlanco, E-28049 Madrid, Spain 
^Physics Department, Theory Division, 
CERN, CH-1211 Geneva 23, Switzerland 

Abstract 

In view of the recent resuhs from the MiniBooNE experiment we revisit the global neutrino 
oscillation fit to short-baseline neutrino data by adding one or two sterile neutrinos with eV-scale 
masses to the three Standard Model neutrinos, and for the first time we consider also the global fit 
with three sterile neutrinos. Four-neutrino oscillations of the (3+1) type have been only marginally 
allowed before the recent MiniBooNE results, and become even more disfavored with the new data 
(at the level of 4cr). In the framework of so-called (3+2) five- neutrino mass schemes we find severe 
tension between appearance and disappearance experiments at the level of more than 3cr, and hence 
no satistfactory fit to the global data is possible in (3+2) schemes. This tension remains also when 
a third sterile neutrino is added, and the quality of the global fit does not improve significantly in 
a (3+3) scheme. It should be noted, however, that in models with more than one sterile neutrino 
the MiniBooNE results are in perfect agreement with the LSND appearance evidence, thanks to 
the possibility of CP violation available in such oscillation schemes. Furthermore, if disappearance 
data are not taken into account (3+2) oscillations provide an excellent fit to the full MiniBooNE 
spectrum including the event excess at low energies. 
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I. INTRODUCTION 



Recently first results from the MiniBooNE (MB) experiment [l|, ^ at Fermilab have been 
released on a search for i/^ v^. appearance with a baseline of 540 m and a mean neutrino 
energy of about 700 MeV. The primary purpose of this experiment is to test the evidence of 



Ve transitions reported by the LSND experiment at Los Alamos [3j with a very similar 
LjE range. Reconciling the LSND signal with the other evidence for neutrino oscillations 
is a long-standing challenge for neutrino phenomenology, since it requires a mass-squared 
difference at the eV^ scale, at odd with the values needed to explain atmospheric 0], long- 
[s], 0], solar 7,0, ol, and long-baseline reactor 10| neutrino data. 



baseline accelerator 



It turns out that introducing a fourth (sterile) neutrino [11 
factory description of all data in terms of neutrino oscillations 

constraints from atmospheric 0], solar [ol, and null- result short-baseline (SBL) experi 



ments 



does not lead to a satis- 
because of tight 
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3| (see Ref. [18| for early four-neutrino analyses considering LSND, 
and Refs. [19|, l20|] for recent updates). So-called (2+2) schemes are ruled out by strong 
constraints on a sterile neutrino component in solar as well as in atmospheric neutrino os- 
cillations 2l|] at high significance. Therefore, we will not consider such schemes in the 



following. Also (3-1-1) schemes suffer from a well-known tension 



3etween the LSND appear 



22, 



24, 



25, 



26,^ 



28 . We 



ance signal and null-result SBL disappearance experiments 
will show that recent MB results further aggravate this tension and hence (3-1-1) schemes 
get even more disfavoured. In Ref. 29| a five-neutrino mass scheme of the (3+2) type has 
been considered, arguing that the disagreement between LSND and null-result experiments 
becomes somewhat relaxed compared to (3+1), see also Ref. 26|]. Here we reconsider this 
possibility in view of the recent MB data. Furthermore, we investigate the impact of adding 
a third sterile neutrino on the quality of the global fit. Since we know that there are three 
active neutrinos, the possibility of three sterile neutrinos is appealing for aesthetical reasons. 



Apart from sterile neutrino oscillations, various more exotic explanations of the 



signal have been proposed, for example, neutrino decay (30|, iMi, CPT violation 



32 



,SND 



vio- 



lation of Lorentz symmetry 33], a lepton number violating muon decay fsJ], CPT- violating 
quantum decoherence 351], mass- varying neutrinos 361], or shortcuts of sterile neutrinos in 
extra dimensions [37]. 

In this work we concentrate on the oscillation framework including one or more sterile 
neutrinos at the eV scale. Such models do have an impact on cosmology 38] . First, the 
1 contribute to the effective number of neutrino species at Big Bang 



sterile neutrinos wi 



nucleosynthesis [23 



39!], and second, these models are subject to strong bounds on the sum 



of the neutrino masses in the sub-eV range from the combination of various cosmological data 
sets (see e.g.^ Ref. 40|] for recent analyses). In order to reconcile such neutrino schemes with 



42 



cosmology some non-standard scenario has to be invoked, see for example Refs. [41j, 

The outline of the paper is as follows. In Sec. [Tllwe consider (3+1) four-neutrino schemes, 
and give some details on the used data and their analysis. Sec. [Till is devoted to (3+2) five- 
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neutrino schemes, discussing the compatibihty of LSND and MB in such schemes, as well 
as the problems of these models to reconcile appearance and disappearance experiments. In 
Sec. lIVI we extend the (3+2) scheme by adding a third sterile neutrino to a (3+3) six-neutrino 
model and investigate whether the global fit improves significantly. We summarise in Sec. |Vl 
In appendix|A]we discuss the mechanism to reconcile LSND and MB by CP violation in (3+2) 
schemes, appendix |B] contains a parameter counting of general sterile neutrino oscillation 
schemes, and in appendix Owe consider details of the analysis of atmospheric neutrino data 
in such models. 



II. (3+1) FOUR-NEUTRINO MASS SCHEMES 

The (3+1) four- neutrino spectra are a small perturbation to the standard three-active 
neutrino case. A cluster of three neutrino mass states accounts for the "solar" (Arrigi) and 
"atmospheric" {Am^^) mass splittings. The fourth mass state is separated by an eV-scale 
mass gap to account for the LSND oscillations, and there is only small mixing of active 
neutrinos with this mass eigenstate. 



A. Appearance data in (3+1) schemes 



We start our discussion of the (3+1) mass scheme by considering the SBL appearance 
experiments in the — » (or z/^ — Ue) channel, including the LSND evidence, the 
bounds from KARMEN 3 and NOMAD ij, and the recent MB data. A combined 
analysis of LSND and KARMEN can be found in Ref. 4J|. In the approximation Am^^ ~ 
Am^i ~ 0, the SBL appearance probability in (3+1) schemes is equivalent to the two- 
neutrino case, where the effective mixing angle is determined by sin^26'sBL = 4|?7e4p|?7^4p. 
Therefore, the analysis performed by the MB collaboration [l|, y] directly applies to (3+1) 
schemes. We comment only briefly on this case, with the main purpose to check our analysis 
against the official MB results. 

For our re-analysis of LSND we fit the observed transition probability (total rate) plus 
11 data points of the L/E spectrum with free normalisation, both derived from the decay- 
at-rest data [3]. For KARMEN the data observed in 9 bins of prompt energy as well as the 
expected background IJ] is used in the fit. Further details of our LSND and KARMEN 
analyses are given in Ref. 3l|- For NOMAD we fit the total rate using the information 
provided in Ref. 15||; our exclusion curve is in good agreement with the result presented in 
that reference. 

The MB re-analysis is based on the z/^ neutrino fiux, efficiencies, and energy resolution 
provided in Ref. [2I, folded with the charged- current quasi-elastic (CCQE) cross section, 
to obtain a prediction for the CCQE event excess from z/^ v^. oscillations. We calibrate our 
simulation to the official MB analysis using the prediction for two example points provided 



3 



10 



^ 10 



10" 



10 



I I I 1 1 llj 1 1 1 11 I 1 1 llj 1 — I I I 1 1 III 1 — I I I 1 1 1 1 



MB475C 



LSND + KARMEN 
NOMAD 




90%, 99% CL 

±1 I I 



_i I 



J I 



J 



lO"* 10"^ 10"^ 10"^ 10° 



sin^20 



SBL 



FIG. 1: Allowed region for MB475 (solid and dashed curves) and LSND+KARMEN+NOMAD (shaded 
regions) at 90% and 99% CL (2 dof) in (3+1) mass schemes. 



in Ref. pU]. For the fit the spectrum of excess events binned in reconstructed neutrino 
energy from Fig. 2 of Ref. [l| is used, where the error bars include statistical errors and 
the uncertainty from the background prediction. Detailed technical information the MB 
oscillation analysis is available at the webpage Ref. j45|, including efficiencies and error 
correlations. Our MB results are in good agreement with the official MB analysis as described 
in that webpage. 

MB data are consistent with zero (no excess) above 475 MeV, whereas below this energy a 
3.6cr excess of 96±17±20 events is observed. Whether this excess comes indeed from z/^ z/g 



transitions or has some other origin is under investigation p|. Lacking any explanation in 
terms of backgrounds or systematical uncertainties we take these data at face value, and 
in some cases we will use all 10 bins of the full energy range from 300 MeV to 3 GeV in 
the fit ("MB300"). However, as discussed in Refs. l2|, two-neutrino oscillations cannot 
account for the event excess at low energies. We confirm that the quality of the (3+1) MB 
fit drastically worsens when the two energy bins between 300 and 475 MeV are included in 
the fit. Therefore, we follow the strategy of the MB collaboration and restrict the (3+1) 
analysis to the energy range from 475 MeV to 3 GeV ("MB475"). 

The bound from MB475 data is shown in Fig. [1] in comparison with the allowed region 
from the combined LSND, KARMEN, NOMAD data. In agreement with Refs. [l|, 13] we find 
that the 90% CL regions do not overlap. A marginal overlap appears if both data sets are 
stretched to the 99% CL. If all data are summed we find a best fit point with Xmin = 26.6 for 
(29 — 2) degrees of freedom (dof). Although this leads to a very good nominal goodness-of-fit 
(gof), the figure clearly shows that there is significant tension between MB and LSND. A 
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powerful tool to evaluate the compatibility of different data sets is the so-called parameter 
goodness-of-fit (PG) criterion discussed in Ref. j46|. It is based on the function 

2 _ 2 _ 2 /I \ 

XPG Xtot,min / j Xi,min ' \ / 

i 

where Xtot,min the minimum of all data sets combined and Xi,min is the minimum of the 
data set i. This function measures the "price" one has to pay by the combination of the 
data sets compared to fitting them independently. It should be evaluated for the number 
of dof corresponding to the number of parameters in common to the data sets, see Ref. [3] 
for a precise definition. Applying this test to check the compatibility of MB with the other 
SBL appearance data, we find Xpg — '^■^ (2 dof), corresponding to a PG of 2.5%. If we 
test the compatibility of the three data sets MB, LSND, and KARMEN+NOMAD we find 
XpQ = 13.7 (4 dof), and PG = 0.8%. In the latter case the slight tension between LSND 
and KARMEN also contributes to the x^? whereas in the first case this tension is removed 
since they are added into one single data set. 



B. Global SBL data in (3+1) schemes 



Now we proceed to the global four-neutrino analysis, addin g al so the information from 
disappearance experiments. We include the Bugey [l7|, Chooz 47|], and Palo Verde 48[ Ue 
reactor experiments, as well as the CDHS |l6| u,, disappearance experiment. Details of our 
Bugey and CDHS fits can be found in Ref. [271]. Furthermore, atmospheric neutrino data 
give an important constraint on the parameter (25I . 28], which is equal to |f//i4p in (3+1) 
schemes, see Ref. 49|] and appendices [HI O This information is crucial for the region of 
Am4i ^ 1 eV^, where the bound from CDHS disappears. We use the updated atmospheric 
neutrino analysis from Ref. 2Q|, which includes also recent K2K |^ and MINOS ^ data, 
and include the bound on rf^ as one single data point in the SBL fit. The total number of 
data points in the (3+1) analysis is 



A^APP — ll(LSND) + 9(KARMEN) + l{NOMAD) + 8(mB475) — 29 
A'dIS = 60(Bugoy) + l(Chooz) + l(Palo Verde) + 15(CDHS) + l(ATM) = 78^ 

It is well known that (3+1) schemes suffer from a tension between the 



tot 



107. 



(2) 



SND appearance 



signal and the bounds from disappearance experiments, see, e.g., Refs. |22l. l23l. I24J. l25l. l26l. l27| . 
Reactor experiments constrain the parameter |f/e4p ^ 1, CDHS and atmospheric neutrinos 
limit |f//x4p -C 1, whereas the LSND oscillation amplitude is given by 4|[/e4p|f^/i4p- Because 
of this tension (3+1) schemes have been already disfavoured before the recent MB results, 
see Ref. [isj. Testing the compatibility of LSND with all the no-evidence experiments (NEV) 
without MB leads to a Xpc = 20.9 for 2 dof, which indicates an inconsistency at high CL. 
Since MB475 data is also in conflict with the LSND signal the new data adds to the tension 
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FIG. 2: Allowed regions in (3+1) schemes from no-evidence (NEV) data including MB475 (solid and 
dashed curves) and LSND (shaded regions) at 90% and 99% CL (2 dof). 



and we find 

XpG = 24.7 (2 dof) , LSND vs NEV (incl. MB475) . (3) 

Alternatively one may test tlie compatibility of the three data sets LSND, 
MB475+KARMEN+NOMAD (NEV-APP), and disappearance experiments, which gives 

XpG = 24.8 (4 dof) , LSND vs NEV-APP vs DIS . (4) 

These numbers formally imply inconsistency at more than 4a, and hence we conclude that 
(3+1) schemes are highly disfavoured by recent data. This level of incompatibility is alre ady 
close to the one of solar and atmospheric neutrino data in (2+2) four-neutrino schemes 
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19 



20| with a Xpg (2+2) ~ ^0.7 (1 dof). Let us mention that the global best fit point of all 
data has Xmm — 100.7 for (107 — 3) dof, which nominally gives an excellent gof. However, in 
the standard gof test the problems in the fit get diluted by the large number of data points, 
and the disagreement of different data sets becomes only visible when they are analysed 



separately and then compared to each other, see Ref. [46|] for a detailed discussion. 

The bound from NEV data in comparison with the LSND region is shown in Fig. [21 
Let us mention that the shape of the NEV bound does hardly change by adding MB data, 
however the of the global best fit point increases significantly (see above). Note also, 
that for our LSND analysis we use only the decay-at-rest data, where the appearance signal 
is most significant.^ If the global LSND data including also the decay-in-fiight events are 



The LSND allowed region of Fig. [2] consists of a connected band, which shows that the fit is dominated by 
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used, the LSND regions shift to shghtly smaller values of sin^ 25sbl and the disagreement 
with NEV gets somewhat less severe, see Ref. VI\ for a discussion of this issue. 

To summarise, (3+1) schemes get further disfavoured by MB data for two reasons. First, 
in these models LSND and MB are in disagreement at the level of 90% to 98.5% CL 
and second, the tension between LSND and NEV data gets worse due to MB. 



III. (3+2) FIVE-NEUTRINO MASS SCHEMES 



Five-neutrino schemes of the (3+2) type are a straight-forward extension of (3+1) 
schemes. In addition to the cluster of the three neutrino mass states accounting for "solar" 
and "atmospheric" mass splittings now two states at the eV scale are added, with a small 
admixture of Ve and to account for the LSND signal. In Ref. 29|] it has been argued that 
in (3+2) schemes the tension between LSND and NEV data becomes significantly relaxed 
compared to the (3+1) case. Here we reconsider this possibility in the light of the new MB 
data. 



A. Appearance data in (3+2) schemes 

First we consider appearance data only (LSND, KARMEN, NOMAD, and MB). In the 
SBL approximation /S.m\^ /S.m\^ ^ 0, the relevant appearance probability is given by 

Py.^ue = ^\Ue4\^\U^,4\'^ Sin2 04i+4|f/e5nt/^5r slu^ 

+ 8 \ Ue4U^4Ue5Uf,5\ siu (/)4i siu (/)5i COs(<^54 - S) , (5) 

with the definitions 

<P^J = , ^ = arg {u:,U,4Ue,U;,) . (6) 

Eq. holds for neutrinos (NOMAD and MB); for anti-neutrinos (LSND and KARMEN) 
one has to replace 6 —6. Since Eq. ([5]) is invariant under the transformation 4 5 
and 6 —6, we can restrict the parameter range to Amf^ > 0, or equivalently Amg^^ > 
Amli, without loss of generality. Note also that the probability Eq. ([5]) depends only on 
the combinations lUe^Uf^^l and \Ue5U^5\, and therefore, the total number of independent 
parameters is 5 if only appearance experiments are considered. 



the total rate and the spectral information available to us is not strong enough to produce disconnected 
regions as obtained from an event-based likelihood analysis 44]. However, the location and size of our 
region are in very good agreement with Ref. \4^ . and moreover, the regions of the parameter space which 
are relevant for the combination with NEV data are reproduced very well. 
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FIG. 3: MB spectral data in bins of reconstructed CCQE neutrino energy. The histograms show the 
prediction at the best fit points in (3+2) mass schemes for SBL appearance data LSND, KARMEN, NOMAD, 
MB (left), and for the global data (right). For the solid histograms the full MB energy range has been used 
in the fit (MB300), whereas for the dashed histogram the two lowest energy data points have been omitted 
(MB475). The corresponding parameter values are given in Tab. IH 

Non-trivial values of the complex phase 6 lead to CP violation, and hence in (3+2) 
schemes much more flexibility is available to accommodate the results of LSND (anti- 
neutrinos) and MB (neutrinos).^ Indeed we find that MB is perfectly compatible with 
LSND in the (3+2) framework. In Fig. [3] (left) we show the prediction for MB at the best 
fit points in the combined MB, LSND, KARMEN, NOMAD analysis. Clearly MB data can 
be fitted very well by simultaneously explaining the LSND evidence; we have checked that 
the prediction for the LSND oscillation probability is within the la range of the observed 
value. In this case also the low energy MB data can be explained, and therefore, in contrast 
to (3+1) schemes, (3+2) oscillations offer an appealing possibility to account for this excess. 
In the following we will present results from both MB data sets, MB475 as well as MB300. 
Note that for MB475 the number of data points used in our analysis is given in Eq. ([2]), 
whereas for the case of MB300 two more data points should be added to iVApp and iVtot- The 
parameter values and the minima at the best fit points are given in Tab. [B In both cases, 
MB475 and MB300, a gof of 85% is obtained, showing that MB is in very good agreement 
with global SBL appearance data including LSND. 

In Fig. m (bottom) the is shown as a function of the CP phase 6. The data prefer 



^ The possibility to use CP violation to reconcile LSND with a possible null-result of MB neutrino data was 
pointed out in Ref. 31 1 in the framework of neutrino decay, and later in Ref. in relation with (3+2) 
oscillations. 
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data set 










5 


Xmin/dof 


gof 


appearance (MB475) 


0.044 


0.66 


0.022 


1.44 


1.127r 


16.9/(29 - 5) 


85% 


appearance (MB300) 


0.31 


0.66 


0.27 


0.76 


l.Olvr 


18.5/(31 - 5) 


85% 


















global data (MB475) 


0.11 0.16 


0.89 


0.12 0.12 


6.49 


1.647r 


94.5/(107 - 7) 


63% 


global data (MB300) 


0.12 0.18 


0.87 


0.11 0.089 


1.91 


1.447r 


104.4/(109 - 7) 


41% 



TABLE I: Parameter values, x^, and gof of the best fit points for SBL appearance data from LSND, 
KARMEN, NOMAD, MB (upper part), and global data (lower part) in (3+2) schemes. Mass-squared 
differences are given in eV^. Results are shown without (MB475), and including (MB300) the low energy 
data from MB. 



values in the range vr < 5 < 27r in order to reconcile LSND and MB. However, as visible in 
the figure no pronounced minimum appears and a rather broad range of 5 values lead to a 
good fit, including also values rather close to the CP conserving value 5 = vr. For MB300 
the best fit even occurs at 5 = l.Olvr, compare Tab. [H In appendix|X]we give an explanation 
how LSND and MB can be reconciled with values of 5 that close to tt. As discussed there, 
this mechanism is based on a delicate interplay of the three terms in the probability of 
Eq. ([5]), and a "tiny amount of CP violation" suffice to make LSND and MB compatible. 
However, this mechanism requires rather large values of the appearance amplitudes \UeAU^A\ 
and \Ue^U^^\, which are not compatible with disappearance data, and hence solutions with 
5 close to vr are not possible in the global fit (see upper panel of Fig. H] and the discussion 
in the following subsection). 

The allowed range in the plane of the two mass-squared differences is shown in Fig. |5l 
Again we observe that the solution is not particularly fine-tuned and a rather wide 90% CL 
region appears. The fit using MB300 is somewhat more constrained by the requirement to fit 
the low energy excess in MB. The fact that for MB300 the best fit occurs for Am^^ Am\i 
is not statistically significant; the allowed region extends far into the hierarchical regime. 

B. Global SBL data in (3-|-2) schemes 

Now we proceed to the global analysis in (3+2) schemes, to see whether the successful 
description of all appearance data found in the previous sub-section can be reconciled also 
with the bounds from disappearance experiments. The (3+2) survival probability in the 
SBL approximation is given by 

P,^^,^ = 1 - 4 I 1 - 5^ |f/,,n ^ |f/„,|2 sin^ 0,1 - 4 |f/„4r|t/a5|' sin^ (7) 

\ j=4,5 / i=4,5 

where (pij is given in Eq. (I6l). Similar as in the (3+1) case, also for (3+2) schemes atmospheric 
neutrino data provide an important constraint on oscillations with sterile neutrinos. The 
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FIG. 4: The in (3+2) schemes as a function of the CP phase 5 defined in Eq. ^ for appearance data 
from LSND, KARMEN, NOMAD, and MB (bottom), and for global data (top). Results are shown without 
(MB475), and including (MB300) the low energy data from MB. All other parameters have been minimised, 
respecting the constraint Am^j^ > Am|j^. 

five-neutrino atmospheric neutrino analysis is discussed in detail in the appendices [B] and [Cl 
It turns out that the same constraint Xatm('^m) four- neutrino case applies, where 

now the definition = lU^^^l"^ + |?7^5p has to be used (see appendix [Bl). 

The results of our global (3+2) fit are summarised in the lower part of Tab. [H where 
the parameter values, the cind the gof of the best fit points are given, again for both 
MB options, MB475 and MB300. The allowed regions for the global fit in the plane of the 
mass-squared differences are shown in Figs. [6] and [71 

In the right panel of Fig. [6] we show the projections for (3+1) and (3+2) schemes. 
Comparing the two best fit points provides a method to assess the relative quality of the 
fit in the two models (likelihood ratio test). We find that introducing the second sterile 
neutrino leads to the relative improvement of the fit of 

Xmin, global (3+1) ~ Xmin, global (3+2) =6.1 (4 dof) , (8) 

where the number of dof corresponds to the additional 4 parameters introduced by moving 
from (3+1) to (3+2). Hence, (3+1) can be rejected only at the 81% CL with respect to 





1 1 1 










1 


1 


1 1 












MB 300 










1 










- 


— 


, — 
















— 1 — 






























1 




























































1 




























































































































































^ 






global data 




























1 1 1 








1 






1 1 1 









1 1 1 1 


































B30C 








f — 










































































































> 










































N 


MB475 / 
























-t 














y 
























/- 
























































































































ippearance 


data 












































1 


1 













10 
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FIG. 5: Allowed regions for SBL appearance data in (3+2) schemes at 90%, 95%, 99%, 99.73% CL (2 dof) 
in the plane of Am|]^ and Amgj^. All other parameters have been minimised. We use data from LSND, 
KARMEN, NOMAD, and MB475 (left) or MB300 (right). 



the (3+2) model. This explains also the "stripes" at Am|]^ and Amg^^ around 1 eV^, which 
appear at 99% CL in Fig. El They correspond to the (3+1) solution, which is always present 
as limiting case in (3+2). Also note, that in the case of appearance data alone we find 

Xmin, APP(3+1) ~ Xmin, APP (3+2) =9.7 (3 dof) . (9) 

Comparing this number with Eq. ([8]) shows that the main improvement in (3+2) schemes 
is to reconcile LSND and MB, whereas it is not possible to evade efficiently the constraints 
from disappearance data. This result is somewhat in disagreement with the conclusion of 
Ref. [is^. A possible source of this different result might be the inclusion of atmospheric 
neutrino data in the fft, which is quite important to constrain sterile oscillations in the 
sector. Our results are in accordance with the arguments presented in the appendix of 
Ref. (26i. 

In the upper panel of Fig. H] the of the global fit is shown as a function of the complex 
phase 5. One can see from that figure that the global data prefer values of 5 close to 
"maximal" CP violation at 5 = 3-71/2. The best fit values are 1.647r and lAAir for MB475 
and MB300, respectively. As discussed in appendix |Al reconciling MB and LSND with 
values of S very close vr (as found for appearance data only) requires rather large values of 
the appearance amplitudes \Ue4Uf^4^\ and \Ue5U^5\, close to the upper bound from unitarity. 
Such large values are not compatible with disappearance data, and hence solutions with 6 
close to TT are not possible in the global fit. 

From the values given in Tab.Uit appears that the model provides a very good fit to the 
data. However, as in the (3+1) case the problem appears when the compatibility of different 
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FIG. 6: Left: Allowed regions for global data in (3+2) schemes at 90%, 95%, 99%, 99.73% CL (2 dof) in 
the plane of Arn^i and Arri^i, with all other parameters minimised. Right: projected onto the Am^i 
axis, using Am|]^ < Am^i. For comparison we show also the projection in (3+1) schemes. The vertical 
dashed lines indicate the of the (3+2) and (3+1) best fit points. For MB the two lowest energy data 
points have been omitted (MB475). 



data sets is considered. Let us divide the global data into appearance and disappearance 
experiments and check their compatibility with the PG test [46| according to Eq. ([T]). We 
find the following Xpc values, for global data without MB, with MB475, and with MB300:^ 

xIg = 17.5 , PG = 1.5 X 10"^ (no MB) 
APP vs DIS: Xpg = 17.2 , PG = 1.8 x 10~^ (MB475) (10) 
Xpg = 25.1 , PG = 4.8 X 10"^ (MB300) 



The PG values have been calculated for 4 dof 46|] . This number corresponds to the four inde 



pendent parameter (combinations) {Ue^U^^l, \ Ue5U^5\, Amli, Amg^, representing the minimal 
number of parameter (combinations) in common to the two data sets. From Eq. flTI]]) we 
conclude that also in (3+2) schemes a severe tension exists between appearance and disap- 
pearance experiments. If MB475 is used the result is very similar to the situation without 
MB data implying inconsistency at about 3.1cr, whereas in case of the full MB data the ten- 
sion becomes significantly worse (about 4(t), since appearance data are more constraining 
because of the need to accommodate LSND as well as the MB excess at low energies. 



We have tested this result explicitly with the "official" MB analysis available at Ref. [45[. Using the MB 
from that source the PG test for appearance and disappearance data gives Xpg = ^"^-^ (MB475) and 
24.6 (MB300), in good agreement with our results displayed in Eq. ([1 
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FIG. 7: Left: Allowed regions for global data in (3+2) schemes at 90%, 95%, 99%, 99.73% CL (2 dof) in 
tlie plane of Arn^i and Arri^i, with all other parameters minimised. Right: projected onto the Am^i 
axis, using Amli < Amgj^. The vertical dashed line indicates the of the best fit point. For MB the two 
lowest energy data points have been included (MB300). 

The tension between appearance and disappearance data is illustrated in Fig. [HI where 
we show the projections of the allowed regions in the plane of the appearance amplitudes 
\Ue4U^i\ and It/gsf^/isl- The opposite trend of the two data sets is clearly visible, especially 
when the low energy excess in MB is included (right panel). Note that an overlap of the 
regions visible in that figure does not prove that there is indeed an overlap of the allowed 
regions in the full parameters space since only a projection is shown. The "common" values 
in the plane shown in the plot might correspond actually to different locations in the space 
of Amli and Amg^. However, if no overlap is visible in that projection at a certain CL there 
is also no overlap at that CL in the full parameter space. 

Comparing the numbers for MB475 and MB300 given in Eq. (fTOl) it becomes obvious 
that the MB low energy excess is a severe problem in the global (3+2) fit, although a very 
good fit can be obtained for appearance data only. This is also apparent from the Xmm 
values given in Tab. [H Adding the two additional MB data points at low energy leads to an 
increase of the best fit of about 10 units from 94.5 to 104.4. Indeed, using the global data 
the MB excess cannot be fitted, as visible in the right panel of Fig. [3l where we show the 
prediction for the MB spectrum at the global best fit point. The reason is that to explain 
the excess relatively large values of |t/e4t^/t4| and \Ue5U^^\ are required (see Fig. [8], right), 
which are inconsistent with disappearance data. 

Before closing this section we give the results of an alternative consistency test for the 
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FIG. 8: Allowed regions at 90% and 99% CL in (3+2) schemes for appearance data (shaded regions) and 
disappearance data (dashed and solid curves) projected onto the plane of |C/e4C^/^4| and \Ue5U^5\. In the left 
panel the two lowest energy data points in MB have been omitted (MB475) , whereas in the right panel the 
full MB energy range has been used in the fit (MB300). 

(3+2) model. Instead of dividing the global data into appearance and disappearance exper- 
iments, we now consider the two data sets LSND and all the remaining NEV data, similar 
as done in Eq. ([3]) and Fig. [2] for the (3+1) schemes. In case of (3+2) this analysis gives 

Xpg = 21.2 (5 dof), PG = 7.8 x 10"^ {3Aa) LSND vs NEV (incl. MB475) . (11) 

Here 5 dof have been used, corresponding to the 5 parameters in common. Without MB we 
find in this case Xpc ~ 19.6, PG = 1.5 x 10^'^ (3.2o"). Hence, the PG test gives a disagreement 
between LSND and the remaining SBL data similar to the disagreement between appearance 
and disappearance data found in Eq. (fTOll . 



IV. (3+3) SIX-NEUTRINO MASS SCHEMES 



Since there are three active neutrinos it seems natural to consider also the case of three 
sterile neutrinos. If all three additional neutrino states have masses in the eV range and 
mixings as relevant for the SBL experiments under consideration, such a model will certainly 
have severe difficulties to accommodate standard cosmology [40|] , and one has to refer to some 
non-standard cosmological scenario 4l|, |42|, |43| . Here we leave this problem aside and focus 



on neutrino oscillation data, investigating how much the fit improves with respect to the 
five-neutrino case. 
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The relevant oscillation probabilities are easily generalised to the (3+3) scheme: 
Pu^^u, = 4 ^ It/eiHf/^iP sin^ii 

i 

+ 8 ^ \UeiUf,iUejUf,j I sin 0^ sin 0^1 cos(0jj - 6ij) , z, j = 4, 5, 6 , (12) 

i,j<i 

with the definitions 

<Pv = . = arg {U:^U,,U,,U;,) . (13) 

Eq. f|T2l) holds for neutrinos, for anti-neutrinos one has to replace 5ij —>■ —Sij. Note that 
only two phases 6ij are independent. The survival probability is given by 

-4^|f/„,n[/,^fsin2 0,^., 2,j=4,5,6. (14) 

i,j<i 

Atmospheric data is included in a similar way as in the previous cases, by using x^idp.), 
where now we define = ^j=4 5 g If^^iP) see appendix [Bl In general, one of the new 
mass-splittings could also fall into the atmospheric range of fewxlO"^ eV^. We have not 
considered such degenerate situations, and we always assume that all Am^^, Am^^-, with 
i,j = 4, 5, 6 are infinite for the atmospheric neutrino analysis. 

We have scanned the of global data in the 11 dimensional parameter space, where a grid 
of 81 X 81 X 81 values for the mass-squared differences has been used, spaced logarithmically 
from 0.1 to 20 eV^. In each point of this grid the remaining 8 parameters have been minimised 
by a standard optimisation routine. 

The results of our search are shown in Fig. [HI where we display the global as a function 
of Aml^. For comparison also the x"^ curves in the case of (3+2) and for MB475 also for 
(3+1) are shown. Note that different from Sec. Illlt here we do not impose any constraint 
on the ordering of the mass-squared differences, and the full s-dimensional space is scanned 
for the (3+s) scheme. This implies that the result is symmetric with respect to the s mass- 
squared differences, and the functions projected onto any of the {Amf^)-a.xes {i > 4) are 
identical. (Just for convenience we label the horizontal axis in Fig. [9] as Amli, since it is 
available in all schemes.) Furthermore, as a consequence of this symmetry there are 2 (3) 
degenerate minima in the five (six) neutrino case, which corresponds to a re-labelling of the 
mass states. Note also, as it must be, the Xmin of (3+(s — 1)) model is the maximal 
value in (3+s), since the (3+(s — 1)) fit is always available as limiting case. 

The global (3+3) best fit points are summarised in Tab.HTl From the table and Fig. [Hlone 
can see that there is only a marginal improvement of the fit by 1.7 units in for MB475 
(3.5 for MB300) with respect to (3+2), to be compared with 4 additional parameters in 
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horizontal dashed lines indicate the corresponding best fit values. 





Am|i 




Amii 




gof 


^(3+2) ^(3+3) 


CL 


MB475 


0.46 


0.83 


1.84 


92.8/(107 - 11) 


57% 


1.7/4 


20% 


MB300 


0.46 


0.83 


1.84 


100.9/(109 - 11) 


40% 


3.5/4 


52% 



TABLE II: Best fit points for global data in the (3+3) scheme. Mass-squared differences are given in 
eV^. We give also the difference between the (3+2) and (3+3) best fits. The last column shows the CL 
at which (3+2) is accepted with respect to (3+3), as derived from the x^ difference evaluated for 4 dof, 
corresponding to four additional parameters in the (3+3) model. 

the model. Hence, we conclude that there are no qualitatively new effects in the (3-1-3) 
scheme. The conffict between appearance and disappearance data remains a problem, and 
the additional freedom introduced by four new parameters does not relax significantly this 
tension. 



V. SUMMARY 

We have considered the global fit to SBL neutrino oscillation data including the recent 
data from the MiniBooNE (MB) experiment [l| in the framework of four-, five-, and six- 
neutrino oscillations. We have divided the global data into various sub-sets and tested their 
consistency within the sterile-neutrino oscillation framework. These results are summarised 
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X'lobai/dof Xpc/dof 
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Xgiobai/dof xlc/dof 


PG 




DIS vs K+N+L 


95 5/96 


14.8/2 


6.1 X 10"^ 


92.1/92 


17.4/4 


1.5 X 10" 


-3 


DIS vs K+N+L+MB475 








94.5/100 


17.2/4 


1.8 X 10" 


-3 


DIS vs K+N+L+MB300 








104.4/102 


25.1/4 


4.8 X 10^ 


-5 


DIS vs K+N+MB475 


70.5/93 


0.1/2 


0.95 


68.9/89 


1.1/4 


0.89 




DIS vs K+N+MB300 








79.1/91 


10.3/4 


3.6 X 10" 


-2 


NEV vs L 


95.5/96 


20.9/2 


2.9 X 10-5 


92.1/92 


19.6/5 


1.5 X 10" 


-3 


NEV+MB475 vs L 


100.7/104 


24.7/2 


4.3 X 10-^ 


94.5/100 


21.2/5 


7.8 X 10- 


-4 



TABLE III: Summary table for various consistency checks within the (3+1) and (3+2) schemes. The 
PG test Q as been defined in Eq. ([1]). In the first column, where we give the data sets tested against each 
other, we use the following abbreviations: K (KARMEN), N (NOMAD), L (LSND), DIS corresponds to the 
disappearance experiments summarised in Eq. ([2]), and NEV=DIS+K+N. Results for the (3+3) scheme are 
qualitatively similar to (3+2). 



in Tab. IIIII for the (3+1) and (3+2) schemes. Clearly, in all cases we find severe tension 
between different sub-samples of the data, with the only exception when LSND and the 
low-energy excess in MB are left out, and in this case indeed no sterile neutrinos are needed 
and the standard three active neutrino scheme (3+0) provides a perfect fit to the data. 
Let us summarise our findings: 

1. (3-|-l) four- neutrino schemes are strongly disfavoured because 

(a) recent MB data is incompatible with LSND at the 98% CL [l|, 

(b) the tension between LSND and NEV SBL data becomes more severe due to MB, 
there is no overlap of the allowed regions for NEV and LSND at 99% CL, and 
the PG test implies inconsistency at the level of 4a, 

(c) it is not possible to account for the low energy event excess in MB. 

2. (3-|-2) five-neutrino schemes 

(a) do provide a good fit to LSND and the recent MB data, 

(b) they can account for the low energy event excess in MB, however 

(c) there is significant tension between appearance and disappearance data (accord- 
ing to the PG test at the level of 3a for MB475 and Aa for MB300). 

3. (3-|-3) six-neutrino schemes do not offer qualitatively new effects, the global ii^" 
proves only by about 1.7 (3.5) units for MB475 (MB300) with respect to (3+2), and 
hence, the conflict between appearance and disappearance data remains. 
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The points 2a and 2b might be considered as an interesting hint in favour of (3+2) schemes. 
Since the combined fit of LSND and MB is based on a non-trivial complex phase which 
introduces a difference in neutrino and anti-neutrino oscillations, these results would repre- 
sent the first indication of CP violation in neutrino oscillations. This hypothesis could be 
tested by MB anti-neutrino data, which is currently being accumulated. However, point 2c 
is a challenge for the (3-1-2) model. The conclusions of 2c and 3 strongly rely on the disap- 
pearance experiments Bugey and CDHS. A crucial check would be the confirmation of v^. or 
disappearance at the 1 eV^ scale. Hence, it might be worth to investigate the possibility 
to obtain such information at future reactor experiments 51|, from disappearance experi- 
ments based on low-energy neutrinos from radio-active sources l52|] , or at the near detector 
complex of up-coming long-baseline superbeam experiments (531]. A characteristic signal 
of sterile neutrino oscillations could be obtained at experiments exploring neutral-current 
detection [54 1. 
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APPENDIX A: RECONCILING LSND AND MB IN (3-^2) SCHEMES 

In this appendix we discuss in some detail how LSND and MB are reconciled in (3+2) 
schemes exploring CP violation in the appearance probability. In particular, it is intriguing 
that a very good fit can be obtained with a complex phase b very close to the CP conserving 
value TT, compare Fig. HI To understand this effect we show in Fig. [10] a zoom into the 
region around 5 = vr, and we display in addition to the also the values obtained for the 
oscillation parameters. 

Let us consider the probability Py^^y^ given in Eq. ([5]). A non-trivial possibility to 
suppress this probability can be obtained by requiring cos (054 — 5) = —1. Then one has 

Pv^^v, = 4(g4 sin 04i - gs sin (p^if , cos(054 - S) = -1 , (Al) 

with the abbreviation qi = \UeiU^i\. Hence, the probability is small for ^ q^ and 054 <^ 1. 
This is precisely the behaviour shown in Fig. (TO) when 5 approaches vr from above, Am54 
becomes small and the q.i approach each other. Writing 5 = vr + e one has cos (054 — 5) ~ 
— 1 + Oi^cpl^, e^), Eq. ( ]A1[) is valid, and the oscillation probability is suppressed in MB. 
Now the question arises whether large enough values for Pp^^p^ can be achieved in order 
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FIG. 10: Fit to appearance data LSND, KARMEN, NOMAD, and MB475 (left) or MB300 (right) in (3+2) 
schemes. We show the (bottom), the values of Am|]^, Am54 (middle), and the values of qi = \UeiUfj_i\, 
i = 4, 5 (top) as a function of the complex phase d defined in Eq. [6l The horizontal dashed lines in the top 
panels correspond to the maximal value allowed by unitarity of the mixing matrix. 



to explain LSND. The difference of anti-neutrino and neutrino probabilities is given by 

Pp^^p, - Pu^,^u, = 16 q^qn sin 04i sin 05i sin ^54 sin e ^ 16 q^qn sin^(05i) ^546 , (A2) 

where in the last step ^54, e <^ 1 has been used. Since ^54 and e are small, the other factors 
have to be as large as possible in order to get a sufficient probability for LSND. Indeed, for 
Am|^ ^ 1 eV^ one has sin^ 05i ^ 1, and also the q^ grow for e — ^> (see Fig. [TOl) . Once 
the maximal values allowed by unitarity, q^ = q^ = 1/2, are reached the LSND probability 
is given roughly by Pu^~,d^ ~ 4e^, where we used Pu^^ue ~ (in order to explain MB) and 
054 ^ e (in order to have cos(054 — 5) ~ —1). Using the experimental value Plsnd = 0.0026 
one finds that a fit should be possible for e ^ 0.025 ~ O.OOSvr, in agreement with our results. 

The similar structure of the left and right panels of Fig. [10] suggests that this mechanism 
works equally well for MB475 and MB300, and fitting the low energy excess in MB does not 
affect these considerations. Obviously, this explanation is not valid for 6 < ir, since the CP 
asymmetry Eq. (lA2p has the wrong sign to reconcile LSND and MB. As visible in Fig. [TOl 
the fit jumps into a quite different solution, which anyway gives a poor x^- Also, the local 
minimum around 5 ~ 7r/2 visible in Fig. H] for MB475 requires a different explanation in 
order to obtain the correct sign of the CP asymmetry for these values of 6. Let us also 
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mention that quite large values of ^4 and gs close to the unitarity bound do appear in the fit 
for 6 ^ Tc, since only appearance experiments are used. Such large values are not possible if 
disappearance experiments are included, which basically require that each of the \Uei\, \U^i\, 
i = 4,5 has to be small. This is one reason for the difficulties in reconciling appearance and 
disappearance data, in close analogy to (3+1). 

APPENDIX B: OSCILLATIONS WITH EXTRA STERILE STATES 

In this appendix we discuss in some detail atmospheric and short-baseline neutrino os- 
cillations involving extra sterile neutrino states. For definiteness, we will focus on (3+3) 
schemes; expressions for (3+2) and (3+1) models can be easily obtained by dropping all 
terms containing a redundant "6" or "5" index. Let us order the flavor eigenstates as 
{ue, z/^, I'r, ^si, ^S2j ^83) and introduce the following parametrisation for the neutrino mix- 
ing: 

U = i?36 i?35 i?34 -R26 R25 -R24 R23 R16 Rl5 Ru Rl3 Rl2 , (Bl) 

where Rij represents a complex rotation by an angle 6ij and a phase ipij in the ij plane, 
while Rij is an ordinary rotation by an angle 6ij. Note that rotations involving only sterile 
states {i.e., Rw with both £, £' > 4) are unphysical, and therefore we have omitted them 
from Eq. flBl|) . For the general case with s sterile states, it is convenient to choose the 
(s + 2) rotations i?i4 and R2j as real, and the remaining ones as complex. The matrix U 
then includes 3(s + 1) angles and (2s + 1) phases. 

A number of simplifying assumptions can be made in the analysis of short-baseline as 
well as atmospheric and long-baseline neutrino experiments. For short-baseline experiments 
one can neglect the solar and atmospheric mass splittings, Am2i = and Amg^^ = 0. In this 
approximation, the mixing angles 9i2, 6'i3 and 6'23 disappear from the relevant probabilities. 
Furthermore, matter effects can be neglected. Since we do not consider neutral current 
interactions in our analysis, the r neutrino is essentially indistinguishable from the sterile 
states, as it participates neither in the production nor in the detection processes. Therefore, 
all the angles also disappear. So for (3+3) models we are left with an effective mixing 
matrix 

U = R26 R25 -R24 R16 R15 Ru for SBL, (B2) 

which contains six angles and two CP phases. In general, under our approximations SBL 
experiments depend on 2s angles and (s — 1) phases. For example, in (3+2) models we have 
four angles (^14, ^15, 6^24, 6^25) and one phase (v^is), and the matrix elements \Uei\, \Ue5\, 
\U^4\, \Uf^5\, aigiU^JJ^iUe^U*^) used in Sec. UTIlare combinations of these five parameters. 

For atmospheric and long-baseline experiments (K2K and MINOS) we neglect the mix- 
ing of z/g with other neutrino states at the LSND mass-squared splittings, justified by the 
constraint from Bugey. This corresponds to setting all the angles 9u with £ > 4 to zero. In 
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this approximation, the complex phase (fu can be dropped. Therefore, in (3+3) models we 
are left with an effective mixing matrix 



U = Rse R35 R34 R26 R25 R2A R23 Ri3 R12 for ATM and LBL, (B3) 

which contains nine angles and four CP phases. As a general rule, in our approximation for 
ATM and LBL experiments the matrix U contains (2s + 3) angles and (s + 1) phases. 

From Eqs. (lB2p and (1B3P it is straightforward to see that for any number of extra sterile 
states, s, atmospheric and short-baseline experiments are connected only through the s 
angles 6^2^ with £ > 4 (or, equivalently, through the parameters lU^el"^ with £ > 4). Note 
that in our convention all the non- vanishing CP phases are "private" to either short-baseline 
{e.g., ipi5 and y^ie) or atmospheric {e.g., ipis, ^934, (^35 and v^se) experiments. 

Let us now focus on the probabilities relevant for the analysis of atmospheric and long- 
baseline experiments. The Hamiltonian in the flavor basis is: 

H = UAU^ + V , (B4) 

where U is given in Eq. ( |B3!) . A = diag(0, Aml^, Am|i, Aml^, Amf^, . . .)/2E, and V = 
±V2 G F di&g{2Ne, 0, 0, Nn, Nn, . . . )/2. It is convenient to define Usm = R23R13R12 and 
U = U t/sM- Then we can write: 

H = UHU^ with H = UsmAUI^ + UWU . (B5) 

In order to further simplify the analysis, let us now assume that all the mass-squared 
differences involving the "heavy" states Uf with i > 4 can be considered as infinite: 
Amfj, Amj^, —* 00 for any i = 1,2,3 and i, i' > 4. In leading order, the matrix H takes the 
effective block-diagonal form: 

where H^^^ is the 3x3 sub-matrix of H corresponding to the first, second and third neutrino 
states, and A*^"*^ = diag(Am4]^, Amli, . . . )/2E is a diagonal s x s matrix (the matter terms 
in this block are negligible in the limit of very large Amj^,). Consequently, the evolution 
matrix is: 

and S = USU^ . (B7) 

We are interested only in the elements See, Se^, S^e and S^^. Taking into account the 
block-diagonal form of S and the relations Uei = Sn and U^i = f/^3 = 0, we obtain: 

See = "^ii , = U*2 Sl2 , = U^2 S21 , 5*^/^ = |?7^2p'S'22 + | | ^6*^"'^ . (B8) 

£>4 




21 



The expressions for the probabihties, Pc 



|5'a/3p, are straightforward. Defining 



we note that \U^2? + Y.i>a Wt^^? 



1 and that Uat. - 



Uai for £ > 4, so that |f//,2p 



Therefore: 
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fj,e 



p(3) 



^1 dpi) PjiJ , 



Pr 



en 



:^-d,)p!S 



p,, = {i-d,)'p^'j + j2\u, 



(B9) 



dp. 



(BIO) 



l>4: 



where we have used the fact that the terms containing a factor e*^"^ oscillate very fast, and 
therefore vanish once the finite energy resolution of the detector is taken into account. In the 
above expression P^^^ is the effective probability derived from the Hamiltonian H^^\ which 
has an ordinary three-neutrino term Hsu (including the usual charged- current interaction 
term of the electron neutrino) plus a matter term arising from the sterile part of WVU: 

/O \ 



H^^^ = HsM±V2G 



|t/„2p U:2Ua3 

Eqs. f lBlOp and (IBlip are valid for any number of extra sterile states. 



a=sterile 



(Bll) 



APPENDIX C: ROBUSTNESS OF THE ATM-^LBL BOUND ON 



As discussed in Refs. 
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281], the contribution of atmospheric neutrino data to the 



disappearance data set plays a crucial role in rejecting sterile neutrino models. In this 
appendix we re-consider the bound on d^ in (3-1-1) schemes, generalise it to the (3-|-s) case, 
and investigate the impact of some of the adopted approximations. 

1. Decoupling electron neutrinos 



Let us begin by considering the simplified case Am|j^ = and ^13 = 0. In this limit. 



the electron neutrino state is completely decoupled, so that Pee = 1 and Pep = P^ 



fie 



0.' 



Eqs. ( ]B10|) and (]B11|) then reduce to an effective two-neutrino form in the fi — t sector: 



AE 



P^^ = {l-dpfP^^^ + Y,\U> 



COS 2^23 sin 26'23 
sin 26'23 cos 26*23 



£>4 



\u, 



a2\ 



a=sterile 



\u. 



a3\ 



(CI) 



(C2) 



^ Up to now this approximation has been always adopted in the Hterature, and in the foUowing Sec. IC 2l we 
are going to relax it for the first time. 
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FIG. 11: Allowed regions for atmospheric and long-baseline data in (3+1) schemes at 90%, 95%, 99%, 
99.73% CL (2 dot), and projections of Ax^ over ds = |J/s4p and — |C/^t4p. In panel (b) the purple line 
("Real") is obtained by forcing the complex phase 1^34 to 0° or 180°, while for the blue line ("Complex") 
(p34 is left free. In panel (c) the purple line is hidden by the blue one. 



1 numerical analysis. 



19|, I21I, |49|, and are 



In the case of only one extra sterile state it is possible to perform a fu _ 
The details of such an analysis have been widely discussed in Refs. 
summarised in Fig. [TTl As can be seen from this figure, atmospheric and long-baseline 
neutrino data strongly prefer a pure two-neutrino Uj. oscillation scenario, disfavouring 

both a sterile neutrino contribution to the main Aml^^ oscillations (parametrised by dg = 
|?7s4p) and a mixing of z/^ with the heavy mass eigenstate (parametrised by d^^ = lUf^^l"^). 
In this work we are mainly interested in the bound on c?^, since the other parameters are 
"private" to atmospheric and long-baseline data and can therefore be marginalised. 

Let us now turn our attention to five-neutrino models. Even in this case it is possible to 
perform a full numerical analysis, presented in Fig. [T2r a). As mentioned in appendix [HI in 
principle atmospheric and long-baseline data constrain separately |f//x4p and |t//i5p, which 
we parametrise in terms of (i^ and = (|f/^4p — |f/^5p)/(i^. However, as can be seen from 
Fig- [TWa) the allowed region is practically independent of Furthermore, comparing this 
figure with Fig. [TlTc) it turns out that the bound on c?^ in four-neutrino and five-neutrino 
models is practically the same. In other words, the extra freedom which we have in 
schemes with respect to Au ones is essentially irrelevant for the constraint on c?^. 

In order to understand this result, let us go back to Eqs. ( ICll) and ( ]C2[) and consider the 
differences between 4z/ and 5z/ models. They arise from two facts: 

• In 4i/ models we have P^^ = (1 — c?^)^ pji^J + d'j^, so that the "scaling" term {1 — d^y and 
the "constant" term ^^>4 \ U^\'^ in Eq. (ICip are related to each other. Conversely, in 



23 




FIG. 12: (a) Allowed regions at 90%, 95%, 99%, 99.73% CL (2 dof) in the {d^,, plane from the analysis 
of ATM+LBL data in (3+2) schemes, (b) & (c) Impact of different approximations on the atmospheric and 
long-baseline bound on d^, see text for details. 

5z/ models P^^ = (1 — PjfJ + + ^^)/2, so that the two terms are independent. 

(2) 

• The expression for P^^ is different in the two models, due to the different contributions 
to the sterile matter term in Eq. (]C2p . Again, in five- neutrino models we have more 
freedom. 

The relevance of these differences is illustrated in Fig. [T2r b). All the lines correspond to four- 
neutrino models. The blue line ("Exact") is the exact bound on from atmospheric and 
long-baseline data, and coincides with the one shown in Fig.[TT](c). The red line ("Vacuum") 
is obtained by neglecting the last term in Eq. (lC2p . i.e. by considering only the vacuum part 
of Finally, the green line ("Scaling") is obtained by further neglecting the constant 

term in Eq. ([CT]), thus setting P^^ = (1 - d^f P£^ As can be seen, sterile-induced matter 
terms are completely irrelevant for what concerns the bound on c?^, and the constant term 
in the expression for P^^ plays only a minor role. The real bound arises from the scaling 
term in P^^, which is the same in four-neutrino and five- neutrino models. This explains why 
the differences between the two schemes are so small. Similarly, the weak dependence on 

in models arises from the constant term in Eq. (ICip : in particular, both Figs. [T^ a) 
and [T^ b) show that as this term is decreased the quality of the fit gets worse. 

In summary, although in principle atmospheric and long-baseline data could constrain 
separately all the terms in models with extra sterile states, in practice they are only 

sensitive to the sum of these terms, d^ = X]^>4 I^^aj^P- Furthermore, the bound on d^ is 
essentially independent on the number of extra sterile states. The validity of this approxi- 
mation is crucial for the analysis of (3-1-3) models, where an exact treatment would be very 
hard to do due to the large number of parameters involved. 
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2. Including electron neutrinos 

In the previous section we have seen that the bound on ci^ reflects the abihty of atmo- 
spheric and long-baseline data to effectively fix the total normalisation of /i-like events. This 
is possible in spite of the large normalisation errors (of the order of 20%) on the atmospheric 
neutrino fluxes and the neutrino-nucleon cross sections, since the accurate measurements 
of e-like neutrino events provided by atmospheric data allows to effectively resolve these 
uncertainties. In other words, what really matters to determine the bound on dfj_ is the rel- 
ative normalisation of /i-like and e-like neutrino events. This opens up the possibility that 
subleading effects modifying the distribution of electron neutrino data may have a sizable 
impact on the bound on d^. In this section we investigate in detail this possibility. 

In the context of sterile neutrino schemes there are two types of contributions which 
affect e-like events: (i) "usual" three-neutrino effects induced by Am^i or by 6'i3, and (ii) 
genuine stehle-u effects induced by non-zero On (with £ > 4). The formalism to study the 
three-neutrino effects has been developed in appendix [Bl Following the results of Sec. IC 11 
we assume that the sterile-induced matter effects in Eq. ( IBlip can be neglected, in which 
case the effective Hamiltonian H'^^^ and the corresponding probabilities Pj^^ reduce to the 
usual three-neutrino expressions. Note that again the relevant probabilities are related to 
the three-neutrino ones only through the parameter d^, and apart from the "constant" term 
in the expression for P^^ in Eq. (IBlOp (see previous section) these formulas are completely 
independent from the number of extra sterile species. The results of our analysis are sum- 
marised in Fig. [T2](c). where we compare the d^ bound for the case when the electron is 
decoupled (red line) with the same bound including also three-neutrino effects due to 6'i3 
and Am2i (black line). Note that the Chooz experiment is also included in the fit. As 
can be seen from this figure, an accurate treatment of subleading three-neutrino effects in- 
deed weakens the bound on d^, however the effect is very small and has no impact on the 
conclusions of this work. 

Let us now study effects induced by non-zero values of On. Following the derivation of 
appendix [Bl it is easy to see that in this case U^i and [/^3 no longer vanish, in which case 
Eq. ( IBIOI) should be replaced with expressions involving not only P^^ but also interference 
terms between different entries of S. However, it is still true that Ue2 = Ues = 0, so that: 

Pee={l-defP^!^ + J2\Uei\\ with d, = J2\Uei\'. (C3) 

£>4 £>4 

Motivated by this result, we approximate the effects of non-zero 6u by introducing an 
independent scaling factor for electron events: 

iVe(4) = (1 - 24) iV,(0) . (C4) 

If de is left free to vary without any constraint, the impact on the bound on d^ is dramatic, 
since in this case electron events can no longer fix the 20% flux and cross section normal- 
isation uncertainties. However, the value of de is strongly bounded by Bugey. We have 
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performed a combined analysis of atmospheric+LBL and SBL data taking into account that 
both data samples depend on and de (using the approximation described above). We find 
that the final result is practically the same as in our standard case, where the dependence 
of the ATM data on de is neglected. To illustrate this we show in Fig. [TW c). light-blue line, 
also the bound on c/^ from ATM data by adding to the overall ^ term ((ie/0.012)^, which 
simulates roughly the constraint from Bugey, neglecting that it actually depends on Amii. 
From this figure it becomes clear that once dg is limited by the data from Bugey its impact 
on atmospheric+LBL data is very small. 

In conclusion, the atmospheric bound on d^ is robust under our approximations. However, 
as clear from the above discussion it depends on the assumptions about uncertainties on 
quantities (like fluxes or cross sections) affecting the ratio of e-like and /z-like atmospheric 
neutrino event normalisations. 
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